Momentum transfer in non-equilibrium steady states 
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When a Brownian object interacts with non-interacting gas particles under non-equilibrium conditions, the 
energy dissipation associated to the Brownian motion causes an additional force on the object as a 'momentum 
transfer deficit'. This principle is demonstrated first by a new NESS model and then applied to several known 
models such as adiabatic piston for which simple explanation has been lacking. 

PACS numbers: 05.40.-a, 05.70.Ln, 05.20.Dd 



In nonequilibrium statistical mechanics the mechanical 
coupling between a system and environments still remains 
poorly understood. In the Langevin description, the frame- 
work of energetics was developed during the last decade (HQ] 
but there are certainly many aspects which cannot be grasped 
by such a level of description. For example, when a Brownian 
object is not symmetric, such as a cone or wedge shape, its 
asymmetric properties are not fully reflected in the linear fric- 
tion constant or tensor, y, of the Langevin equation because y 
is non-polar. 

Related to this limitation, or due to our lack of comprehen- 
sion about nonequilibrium Brownian motion, there are a class 
of nonequilibrium phenomena which have refused to be un- 
derstood at a fundamental level. An interesting example is 
the adiabatic piston separating two gases of different temper- 
atures under pressure equilibrium |3j-|5]. The laws of ther- 
modynamics cannot tell whether the piston moves or not . 
Feynman 14J] pointed out that the fluctuations of piston's ve- 
locity should be taken into account. However, the Langevin 
description with linear friction falsely predicts zero mean ve- 
locity. The adiabatic piston is, therefore, still listed among 
major unsolved thermodynamics problems @]. This difficulty 
is also shared by some models of Brownian ratchet motors 
working between ideal gas reservoirs J3]. 

A common solution to these problems is to resort to full and 
general microscopic descriptions, such as the molecular dy- 
namic (MD) simulation or master-Boltzmann equation under 
pertinent perturbative approximations @]. These methods are 
effective in predicting the outcome. For the adiabatic piston, 
the MD simulations |@] and perturbative master-Boltzmann 
equation IT ] give quite consistent results showing that 
the piston moves towards the hotter reservoir. For the ratchet 
models the agreement between MD simulation and perturba- 
tive theory is excellent 0]. When higher order terms are taken 
into account, the perturbative theories can tell the effect of the 
shape of Brownian object ||7|| or of the inelasticity of colli- 
sions, called inelastic piston Il2l 1 1311 and their combinations 
[l5ll .Yet. we still have no physical explanation why the 
nonequilibrium processes give rise to a force and what deter- 
mines its direction. 

In this paper, we will develop a general theoretical frame- 
work to answer to this fundamental problem. The key is to 



explicitly take into account the momentum and mass balances 
under nonequilibrium condition, in addition to the energy bal- 
ance considered by the stochastic energetics J2l ■ Briefly, the 
nonequilibrium energy transfer, or dissipation, leads to a de- 
ficiency in the momentum transfer from the environment to 
the Brownian object, while the gas particle (mass) flux is un- 
changed by the dissipation. We shall call this deficiency the 
momentum transfer deficit due to dissipation or MDD, for 
short. We will show that this MDD is expressed as a form of 
nonequilibrium boundary condition for the momentum flow 
[Eq. ([T} below]. With this condition, many nonequilibrium 
problems which have been hitherto solved case-by-case can 
be explained in a unified manner sometimes even at semi- 
quantitative level. 

In the following, we first describe the basic principle. To 
demonstrate the principle we introduce a simple model of 
nonequilibrium steady state (NESS) and its solution. Then, 
we will apply the basic principle to unexplained problems 
such as adiabatic piston in order to show the universality of 
underlying physics. We extend our principle to include weak 
inelasticity of collisions between the Brownian object and gas 
particles. 

In the elementary setup Fig. [Ha), an ideal gas of temper- 
ature T and pressure p fills to the left of the wall. The wall 
is a Brownian object and its velocity fluctuates. However, it 
macroscopically remains at rest. The collisions of the gas par- 
ticles with the wall strictly satisfy the momentum conservation 
but can be either elastic or inelastic. We assume that the en- 
ergy transfer by individual collision is very small so that the 
double-collision by the wall with the same gas particle is neg- 
ligible. More specifically, the mass of the wall, M, and that of 
gas particles, m, are assumed to satisfy e 2 = m/M «: 1. Our 
interest is the force exerted on the wall by the gas or equiv- 
alently the momentum transfer from the gas to the wall. We 
separate this momentum transfer into two parts; one due to 
the incoming particles toward the wall, £ mv m , and the other 
to the outgoing particles from the wall, 2( _OTV out), where the 
sum is taken over a unit time. The sum of the two momentum 
fluxes gives the force on the wall. 

When the wall's microscopic fluctuations are thermally in 
equilibrium with the gas [Fig.[TJa)(left)], the detailed balance 
condition tells that two momentum fluxes should be equal on 
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FIG. 1. (a) When the piston's mean velocity V is zero, the net mo- 
mentum transfer by the incoming particles, 2 mv in , and by the out- 
going particles, 2(-mv 0lll ), per unit time sum up to give the force 
on the piston. (Left) Equilibrium state where no net energy is trans- 
ferred to the piston. (Right) Non-equilibrium case where the energy 
is dissipated at a rate jSL When Jv > 0[< 0], additional force 
Fmdd < 0[> 0] is exerted on the piston, (b) Cooled/warmed Brow- 
nian piston : The piston (thick bar) and object (diamond) are tightly 
bound and are held by a spring. The gas environments have temper- 
atures T and T' and pressures p and p', respectively. 



the time average, and the sum of the two is the hydrostatic 
pressure p times the surface area, L. Note that unlike a simple 
kinetic theory used in elementary textbooks the individual col- 
lisions can transfer energy between the gas and the wall at the 
microscopic level since the wall fluctuates. It is the detailed 
balance that makes the two momentum fluxes identical on the 
average. On the other hand, when the dissipation carries away 
a part of kinetic energy of the gas upon collision to the outside 
of the system at the rate per unit time [Fig. [TJa)(right)], 
the speed of the outgoing particles is, on the average, less than 
that of the incoming ones. Therefore, the momentum transfer 
by outgoing flux, Yj(-mv out ), should be less than the incoming 
one, 2 mvi n , which should not be influenced by the dissipation 
as long as the double collisions are negligible. This reduction 
in momentum transfer is the MDD, and the resulting force on 
the wall is less than that in the equilibrium by the MDD. This 
additional force due to MDD is exactly the point where the 
Langevin equation with linear friction fails to grasp the left- 
right asymmetry of the system. 

To make this principle more concrete and quantitative, we 
first assume elastic collisions between gas particles and the 
wall. We take the thermal velocity v t h = y/k B T /m as a typical 
normal component of the velocity of incoming particles vj„ up 
to a numerical factor (see below). The first part of momentum 
transfer is £ 'Win — '«^th w co i, where a> co \ ^ pLv t h/2 is the col- 
lision frequency on the wall. We denote by v'(< 0) the typical 
normal component of the outgoing velocities v ollt . The second 
part of momentum transfer is then J](-mv out ) ^ m\v'\a> co i. The 
conservation of mass fluxes imposed the common frequency, 
o) co \ for both incoming and outgoing fluxes. 



Now v' is related to v± though the energy balance condi- 

tion, j Vth 2 - f v' — . Here, we assumed that the parallel 
component of the velocity does not contribute to the energy 
loss. Noting |v'| - v t h for weak energy transfer, the left hand 
side can be approximated by v t h('«Vth - m|v'|). Then the MDD 

par unit time is (mv t h - m\v \)o) co i - -^ a , and the net force on 
the wall is 



F - F ecj + Fmdd, 
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where F eq = pL is equilibrium hydrostatic force and the nu- 
merical constant c is 1 in the above semi-quantitative deriva- 
tion. From the view of the gas, Eq. (Q~|i can be considered as 
a boundary condition for the momentum flux. This additional 
force Fmdd induced by dissipation is the main result of the 
present Letter. 

An interesting realization of MDD, which is also a new 
model of NESS, is illustrated in Fig.[TJb). In two dimension, 
a piston with smooth vertical wall as Brownian object is in 
contact with a gas of temperature T and pressure p. Its hor- 
izontal motion is tightly coupled to another object (rhombus) 
immersed in a different gas environment of temperature T'. 
We can show that when the horizontal diagonal 6\\ and vertical 
diagonal t x of the rhombus are made indefinitely large keep- 
ing i\ = 2(\l(\\ constant, the collisional forces from the sec- 
ond gas converges to the ordinary force of Langevin equation, 
that is the frictional, — y'V and random force, ^2y'k B 7 "<$>,, 
with the friction constant y' = -yV/8 p€imv± (p = p/knT) 
and Gaussian white noise 0, with (©/©;-) = 5(t - t') 11611 . 
Therefore, for T < T the thermal contact dissipates energy to 
the second gas without net transport of momentum between 
the two gases on the time average. 

For weak dissipation, the dissipation rate 7l depends on 
the coupling with the environments only through the friction 
constants, y with the first gas and the aforementioned y' with 
the second gas: 
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For later use, we here show a heuristic derivation of Eq. 
Assuming a kinetic temperature of Brownian motion, k B T^ n , 
we can construct dimensionally /*? by the time constant M/y 
and the temperature gap, k B T -k B T^ n , as /*? = (y/M)(k B T- 
kBTfo n ). Applying the same argument for the second bath, 
-(y'/M)(k B T' - fcB7\ in ), and eliminating T^ n , we 
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obtain Eq. For the standard derivation, see Ref. |0]. 

The linear friction, y = c'pLmv±, can be also obtained heuris- 
tically by a Doppler shift of the momentum transfer, where 
c' — y/8/n from the standard gas kinetics. 

Now substituting Eq. (|2]) to (Q]) we have a concrete form 
of MDD and the force in nonequilibrium. We can show that 
the microscopic approach with the master-Boltzmann equa- 
tion gives exactly the same result if we choose c = yfnJE. 
When the wall is 'cooled', i.e. for T < T, the mean posi- 
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tion of the wall in NESS is displaced leftwards relative to its 
equilibrium position , and vice versa. 

The principle (f2]i is also applicable to the case where the 
collision is weakly inelastic. In Fig.QJb) we remove the rhom- 
bus and the second gas environment, and instead assume the 
restitution coefficient e (1 - e <K 1) for the collision between 
the gas particles and the vertical wall. In this case, the dissi- 
pation rate consists of two parts: 
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The 'house-keeping' heat generation lfl7tl . jty hk , is due to in- 
elasticity of individual collisions. The 'excess' dissipation, 



(e) 



, is due intrinsically to the nonequilibrium Brownian 

diss,ex' J ~ 

motion of the wall. If the wall were rigidly fixed, only hk 



is nonzero. In this case the dissipation per collision is mvi/2- 
mv' 2 /2 = (1 - e 2 )mvl/2 and tt = (1 - e 2 )mv 2 J2 X w col . 



^ e 2 )mvll2 W aj^ = {\-e 2 )mvi, 
Noting (1-e) 2 - 2( 1 -e), the same argument leading to Eq. ((TJ 
gives 
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where Fmdd,!* is force due to the 'house-keeping' MDD 
which reduces the force even for a fixed wall. (A sand bag 
will receive less impact than a hard wall by a bullet.) 

The excess dissipation is expressed in terms of the 
aforementioned kinetic temperature k B T kin as J^L — = 
(M/y)(k B T - IcbT^). Upon a binary inelastic collision, the 
velocity of a Brownian object changes in the same way as that 
of elastic collision if the effective mass M e g = 2M/(l + e) 
is used. The Brownian object then obeys approximately the 
Maxwell distribution oc e -M ca v 2 /(2k B T)^ It i mp ii es k B T^ = 
k B T x (1 + e)/2. Therefore, 
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where the friction constant y is the same as before in the low- 
est order in (1 - e). With the numerical factor c - ynJE, we 
recover the microscopic result. When the dominant house- 
keeping MDD is canceled by the same MDD from the other 
sides, as for an inelastic triangular Brownian object 11411 . it 
is the excess MDD that explains the origin of nonequilibrium 
force. 

We have shown that our simple calculation gives the identi- 
cal result as microscopic approaches up to a numerical factor 
of order one. We note that the microscopic approach is still 
needed to find the correct numerical factor for the Brownian 
object of complicated geometry and to find higher order cor- 
rections to the perturbation. However, our main goal is rather 
to show that the principle ([TJ is a general theory of the force 
under nonequilibrium condition. Below we will apply the ba- 
sic schema Fig.[TJto various known cases and show how our 
principle based on the conserved quantities is fundamental to 
understand the phenomena. 



Adiabatic piston ( with elastic wall) fll^l- Hlll : We apply the 
boundary conditions ([TJ) to the both sides of the piston shown 
in Fig |2ja), with an appropriate sign of the forces and dissipa- 
tion rates as well as taking account of different temperatures. 
By the isobaric condition, the equilibrium force F eq on both 
sides cancels. On the side of temperature T the force, F' MDD , 
should contain the dissipation rate /S' = to assure the 

energy conservation. Both Fmdd and F' MDD are oriented to- 
ward hot side (leftward if T > T'), and the total momentum 
balance about the piston is recovered by the frictional force, 
Fmdd + -Fmdd _ (T + YW - ®> wnere V is the steady state 
velocity of the piston. Combining with Eqs. (QJ and (O the 
steady state velocity reads 



V = - 



i | i \ k B T-k B r 

y + y'\ vth v[ h j M(y~ l + y'~ l ) ' 



(6) 



which is identical to the result obtained from the perturbative 
calculation in ifioll with c = (p - pk B T is assumed 

to verify it.) The correction to dissipation due to the 'meso- 
scopic loss' (y + y')V 2 is of higher order by e 2 and, therefore, 
negligible. This remark applies to all other examples below. 

Inelastic piston lfl2[ 131 : A piston of two inelastic faces 
shown in Fig. [2jb) is in a gas of temperature T and pres- 
sure p, and the faces have coefficients of restitution e (left 
face) and e' (right face), respectively, with 1 - e <K 1 and 
1-e' «: 1. The dissipation rate , MDD and /^mdd, on 
each face satisfy Eq. ([TJ. But the MDD on the two faces has 
different signs, and thus the net force arises only when e + e', 
as fMDD| left + ^MDD| right = (e-e')pL/2. By balancing with the 
frictional force, the stationary velocity V to the lowest order 
in 1 - e and 1-e' (therefore |e - e'| <k 1) and in e is 



V = 
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y + y' 2 



-pL 
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where y y' — c'pLmv t h with c' a constant. The result (0 
agrees with the perturbative results Il2[ 1 1 3fl with c' = y/W/n. 
This elementary example shows, however, that our principal 
formula ([TJ is universal whether or not the origin of dissipa- 
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(b) Inelastic piston M 
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FIG. 2. (a) A microscopic "adiabatic" piston of mass M (vertical bar) 
separates two semi-infinite gases of point-like particles with mass 
m{<K M), The two gases have the same pressure, p, but different 
temperatures, T and T . (b) A macroscopic inelastic piston with 
restitution coefficients, e (left surface) and e' (right surface) is in a 
gas. 
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FIG. 3. A Brownian ratchet proposed in 0](top) consists of two 
triangles (total mass M) that translates along the horizontal axis. This 
model can be mapped to the adiabatic piston (bottom) in the limit, 
6>->0. 



tion is kinematical or dynamical, because the momentum con- 
servation is universally valid. 

Ratchet model in two gas environments: Van den Broeck 
et al. J3l proposed and analyzed a series of Brownian ratchet 
models that move horizontally in contact with two ideal gas 
environments at different temperatures T and T . One typical 
example is shown in Fig. [3f top), where we assumed isobaric 
condition only to simplify the argument without loosing the 
essential point. Microscopic methods concluded that it moves 
steadily with the base of the triangle in hotter environment be- 
ing ahead, i.e. leftwards if T > T . Based on our principle, 
the origin of nonequilibrium force is essentially the same as 
the aforementioned adiabatic piston. Intuitively, if we look at 
only the bases of triangles, it already appears identical to the 
adiabatic piston, Fig.[3bottom). In fact, the sides of the trian- 
gle receive more frequent collisions than on the base but with 
much less impact on the horizontal motion. We can rigorously 
show that in the limit of 6 — > (see Fig. [3}, the momentum 
transfer rate on the sides converges to the equilibrium force, 
pL, without fluctuation or frictional velocity dependence [16]. 
Therefore, in this limit, the effect of side surface vanishes and 
the same principle as the adiabatic piston determines the mo- 
tion of the ratchet model. The result agrees with their pertur- 
bative calculation 17[]. 

Inelastic triangle: Costantini et al. Hill studied a variant 
of above ratchet model using a single triangle but with in- 
elastic surface of restitution constant e. In this case, the net 
house-keeping component vanishes, as if the triangle is in a 
hydrostatic pressure, (1 + e)p/2. On the other hand, the ex- 
cess dissipation /V? ex on the side surfaces are less important 
than that on the base, in the way that the contribution by the 
side surfaces vanishes in the limit S — > 0. In this limit, the 
force balance with frictional drag -yV and Eq. (0 yields 
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This result is identical to the one obtained by microscopic ap- 
proach [Till to the lowest order in 1 -<?, if we choose c = y/n/8. 



In summary we have introduced a unified theory on the gen- 
eration of nonequilibrium force as momentum transfer deficit 
due to dissipation. This principle is applied to a new model 
of NESS, named, cooled/warmed piston, as well as to many 
existing models such as adiabatic piston in a unified man- 
ner. What we clarified here is that, while the energetics at 
Langevin level J2l is enough to treat the dissipation, the dissi- 
pation attributed to Brownian motion plays a decisive role fH 
in the force generation through the MDD. As perspectives, the 
MDD should be taken incorporated in the hydrodynamic de- 
scription of adiabatic piston 11811 . It is of interest to generalize 
the present results to interacting gas particles, for examples 
the boundary thermostats lfl9tl as well as to the contact value 
theorem 12011 under nonequilibrium. 
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